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Abstract. After a dynamically active emergence phase, magnetic flux at the solar surface soon ceases to show 
strong signs of the subsurface dynamics of its parent magnetic structure. This indicates that some kind of discon- 
nection of the emerged flux from its roots in the deep convection zone should take place. We propose a mechanism 
for the dynamical disconnection of the surface flux based upon the buoyant upflow of plasma along the field lines. 
Such flows arise in the upper part of a rising flux loop during the final phases of its buoyant ascent towards 
the surface. The combination of the pressure buildup by the upflow and the cooling of the upper layers of an 
emerged flux tube by radiative losses at the surface lead to a progressive weakening of the magnetic field in several 
Mm depth. When the field strength has become sufficiently low, convective motions and the fluting instability 
disrupt the flux tube into thin, passively advected flux fragments, thus providing a dynamical disconnection of the 
emerged part from its roots. We substantiate this scenario by considering the quasi-static evolution of a sunspot 
model under the effects of radiative cooling, convective energy transport, and pressure buildup by a prescribed 
inflow at the bottom of the model. For inflow speeds in the range shown by simulations of thin flux tubes, we find 
that the disconnection takes place in a depth between 2 and 6 Mm for disconnection times up to 3 days. 
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1. Introduction 



Until the advent of local helioseismology, very little was 
known about the subsurface structure of sunspots and ac- 
tive regions. Differences in the measured rotation rates of 
the nonmagnetic plasma and the magnetic 'tracers' were 
taken as evidence that the latter are somehow 'anchored' 
in deeper layers, but no consensus could be reached as 
to where and how this anchoring should precisely take 



place (e.g. 
ll997tlBecl 



D'Silva fc Howardlll994HJavaraiah fc Go khalc 
2000HSchusslejll984HSchusslerlll987l) . The re- 



cent results from local helioseismology indicate significant 
changes of the thermodynamic and/or magnetic proper- 
ties as well as of the flow field at a depth of l ess than 10 



Mm below sunspots jKosovichev et alJ feoOO; 
200lUCouvidat et al.l20ol ' and active regions 



Zhao et al 



Basu et al 



2004HHindman et al.ll2004|) . These findings are consistent 



with the model of a sunspot as a shallow object, which 
breaks apart into a large number of flux tubes not far 
below the surface l|ParkeJll979T: ISpruitJll98ll IChoudhuril 
Il992lh rather than that of a monolithic plug of magnetic 
flux extending deep into the convection zone (see also the 
detai led discussion of sunspot models bv lThomas fc Weissl 
119921) . 



From the theoretical side, analytical studies and nu- 
merical simulations of magnetic flux tubes have provided 
a consistent picture of the formation of active regions 
and sunspot groups, comprising the storage and ampli- 
fication of magnetic flux near the bottom of the convec- 
tion zone as well as the formation of unstable loops and 
their rise through most of the convection zone. Given a 
strongly super-equipartition initial field of the order of 
10 5 G (10 Tesla) at the bottom of the convection zone, 
large-scale properties of newly-emerged active regions and 
sunspot groups, like low emergence latitudes, systematic 
tilt angles (Joy's law), and asymmetric proper motions 
of the two polarities, can be quantitatively explained (see 
revie ws bv lMoreno-Insertidll997l iFisher et alll2000t iFanl 
2004, and further references therein). 

Numerical studies following the evolution of an unsta- 
ble flux loop from its origin at the bottom of the con- 
vection zone only reach until about 10 Mm from the 
surface, where the thin- flux-tube approximation breaks 
down. Realistic simulations of the subsequent phases are 
still too demanding in view of the computing power avail- 
able today. Therefore, the actual emergence of flux in the 
photosphere and the evolution thereafter has not been cov- 
ered by such simulations (except for rather idealized situa- 



2 



M. Schiissler and M. Rempel: Disconnection of sunspots 



tions aiming at descri bing; the effects on the coronal mag- 
netic structure e.g. fran fe Gibson] |2004 iMagaralEooi 
lArchontis et alj[2004jh On the other hand, a wealth of ob- 
servational evidence indicates a remarkable change in the 
dynamical properties of sunspots and active regions from 
' active' to ' passive' evolution shortly after their emergence 
l|Schiisslerl Il987t ISchriiver fc Title! fl999^ . The emerging 
flux initially displays the clear signature of a rising, frag- 
mented flux tube and evo lves according to its internal 
large -scale dynamics (e.g., iMcIntosbl Il98li IStrous et alJ 
Il996l) . However, the expansion of the bipolar region and 
the proper motion of the sunspots with respect to the 
surrounding plasma decays within a few days after emer- 
gence and larger magnetic structures start to fragment 
into small-scale flux bundles, which are largely dominated 
by the local near-surface flows (granulation, supergran- 
ulation, differential rotation, meridional circulation). The 
magnetic flux is then passively transported by these veloc- 
ity fields and becomes dispersed over wide areas 1 , a pro- 
cess well repr esented by the so- c alled surface transport 
models (e.g.. [Wang et ail Il989t Ivan Balleeooiien et alJ 
Il998l ISriiriiverll200lt iBaumann et al.ll2004lb This behav- 
ior is in striking contrast to what would be expected if 
the emerged flux would follow the evolution of its mag- 
netic roots at the bottom of the convection zone: 

Longitudinal drift: The size range of bipolar magnetic 
regions corresponds to azimuthal wavenumbers between 
rn = 10 and 60, while the magnetic instabilities lead- 
ing to loop formation favor va l ues of m = 1 or 2 
llSpruit fe van Ballegooiienl Il982t ISchiissler et"afl 11994 
iFerriz-Mas fc Schiisslerlll99E ). Consequently, the two po- 
larities of an active region should move apart in longi- 
tudinal direction much further than actually observed. 
Moreover, there is no static equilibrium of the subsurface 
vertical 'legs' of an emerged bipolar region that is con- 
nected to an azi muthal field of 10 5 G at t he bottom of the 
convection zone l|van Ballegooiienl Il982h . The horizontal 
component of the magnetic tension force resulting from 
the bends of the flux tube where it turns horizontal sus- 
tains the drift of the two legs (and their associated poles 
at the surface) in opposite azimuthal directions. Thin-tube 
simulations in fact clearly show that the two poles drift 
apart unrestrained and even move around the whole cir- 
cumference of the Sun within a few months. Clearly, no 
such systematic motion is observed. 

Poleward drift: The emergence of a loop cuts off the az- 
imuthal flow nec cessarv for mechanical equil ibrium of the 
submerged field I Moreno-Insertis et alJ Fl9921. As a result, 



1 When the average flux density exceed s a limit of ab out 
100 G, the relatively stable 'plage state ' JSchriiveJlT987l) IS 
maintained as an intermediate stage: the convective pattern is 
severly disturbed by the magnetic field and supergranulation 
is largely absent, so that the dispersal of the flux is temporar- 
ily diminished. Sunspots actively shape the surrounding flow 
structure by developing the outward moat flow, which probably 
temporarily suppresses the flux transport by supergranulation 
and meridional flow. 



the magnetic structure starts to drift poleward in response 
to the unbalanced latitudinal component of the magnetic 
tension force. No such systematic poleward drift is ob- 
served in bipolar regions or sunspot groups. 

Tilt angle: The tilt angle of sunspot groups with respect 
to the East- West direction arises through the action of the 
Coriolis fo rce on a rising and horizont ally expanding flux 
loop (e.g.. iD'Silva fc Choudhurll^.^ . After emergence, 
the expansion stops and the flux tube should relax back to 
its original East- Wes t orientation du ring the lifetime of a 
large bipola r region jF an et al.|| l994lb No such relaxation 
is observed l|T6th fc GerlefeOO^K 

The observed change of the dynamics of magnetic 
structures from 'active' to 'passive' and the points listed 
above clearly indicate that bipolar magnetic regions some- 
how become disconnected from their deeper roots within 
a few days after emergence. Moreover, since the expected 
motions of sunspots would otherwise be quite marked, the 
complete absence of such behaviour indicates that the dis- 
connection process must reliably work, virtually without 
exception. This represents an important constraint. 

To our knowledge, two mechanisms for the discon- 
nection of emerged fl ux have been proposed so far. 
ISchriiver fc Title! ( 1999J) have suggested subsurface recon- 
nection of the two opposite polarities leading to discon- 
nected U- loops in the upper layers. It does not seem ob- 
vious how this mechanism can achieve the required per- 
fect reliability in view of the strong tendency of the poles 
to mov e apart in the lon gitudinal direction (see Sect.F2J). 
Earlier. [Fan et all l(l994l) had given arguments for the ne- 
cessity of disconnection and proposed a mechanism for the 
'dynamical disconnection' of emerged flux by a local loss of 
lateral total pressure balance following the establishment 
of an isentropic hydrostatic equilibrium along magnetic 
field lines. A s discu ssed in Sec. O the process sketched 
by iFan et aT] l)l994|) does not seem to operate sufficiently 
rapidly to provide disconnection within the first few days 
after emergence. Here we propose a variation of the dy- 
namical disconnection scenario, which is based upon the 
strong, buoyancy-driven upflow associated with the final 
phase of the rise of a flux loop through the convection 
zone. We discuss the various scenarios for disconnection 
in Sect. [3 and provide an illustrative quantitative elabo- 
ration of our model in Sect. 03 The results are discussed 
and put into a broader perspective in Sect.^J 

2. Scenarios for the disconnection of emerged flux 

ISchriiver fc Title! l(l999h have suggested that active re- 
gions literally become disconnected from their roots by a 
subsurface reconnection between the opposite polarities in 
their (fragmented) legs; this would create shallow O- or U- 
loops, which then could float freely with the near-surface 
velocity fields. The problem with this concept is that, as 
explained above, the two legs of a bipolar region tend to 
rapidly move apart from each other owing to the geom- 
etry of the rising loop and as a result of the unbalanced 
tension force in the very deep part. Thus is would require 
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a rather organized and converging flow to bring the oppo- 
site fluxes together to reconnect. Furthermore, the more 
separated the two legs already are, the more difficult is it 
to bring the opposite-polarity flux together. Consequently, 
once reconnection has failed to occur in time for whatever 
reason and the two poles continue to drift further apart, 
they cannot be disconnected any more and separate indef- 
initely. To our knowledge, not a single example of such a 
behaviour of a larger active region has been reported by 
ob servers. 

iFan et all 1^)94) have suggested that the establish- 
ment of hydrostatic equilibrium along a flux tube could 
lead to a dynamical disconnection of the emerged part 
from its deeper magnetic roots, which does not explicitely 
require a change of the magnetic topology (i.e., reconnec- 
tion) : as the subsurface part of the emerged flux tube ap- 
proaches hydrostatic equilibrium, it eventually loses lat- 
eral pressure equilibrium at a certain height. This leads 
to a region of weak, passive magnetic field that effectively 
decouples the parts of the flux tube above and below this 
height. 

This mechanism is closely rela ted to the 'explosion' 
of mag netic flux tubes discussed bv lMoreno-Insertis~et" ID 
(f995): a slowly rising flux loop that maintains approx- 
imate hydrostatic equilibrium along the field lines ex- 
periences a sudden catastrophic weakening of the field 
strength at its apex when it reaches a critical height. 
This explosion height depends on the field strength in 
the part of the flux tube that remains at the bottom 
of the convection zone. The plasma is nearly isentropic 
within the adiabatically rising flux loop whereas the en- 
tropy decreases with height in the surrounding supera- 
diabatically stratified convection zone. As a consequence 
of the higher entropy within the magnetic flux tube, the 
internal pressure decreases more slowly with height than 
the external pressure and thus eventually both become 
equal at a certain height above the b ottom of the convec- 
tion zone l|Schiissler fc Rempel20 02) . When the loop apex 
approaches this 'explosion height', it expands drastically 
since pressure balance forces the magnetic field to become 
very small, so that it is n o longer dynamically relevant 
ijRempel fc Schiisslerl l200lT) . From mixing- length models 
of the solar convection zone one finds that a hydrostatic 
flux tube with a field strength of O(10 5 ) G at the base of 
the convection zone should explode at a depth of less than 
10 Mm below the surface. On the other hand, dynamical 
simulations show that rising flux loops with such initial 
field strengths actually traverse this height range without 
exploding because a flow along the field lines keeps them 
sufficiently far awa y from hydrostatic equilibrium (e.g., 
ICa.ligariet a,l.lll99^. 

The problem with the scenario of IFan et alJ l|1994r) is 
the timescale for the establishment of the 'global' hydro- 
static equilibrium required for the disconnection. The as- 
sociated time scale if of the order of the Alfven travel 
time along the circumference of the flux tube: taking an 
Alfven speed of lkm-s -1 (corresponding to 10 5 G at the 
bottom of the convection zone), this corresponds to about 



a month, much too long to avoid an unrealistic evolution 
of the emerged part. 

However, a more detailed consideration of the develop- 
ment in the upper parts of an emerged flux tube reveals 
a much faster and more local route towards dynamical 
disconnection, avoiding the requirement of a global hy- 
drostatic equilibrium. Consider the following situation. A 
rising flux loop has emerged at the surface to form a bipo- 
lar magnetic region. The subsurface part of the loop has 
not exploded at its formal explosion height because the 
strongly buoyant and rapidly rising plasma in the loop 
maintains a super-hydrostatic gradient of the gas pressure 
(i.e., a steeper upward decline of the pressure than in the 
hydrostatic case) co nnected with an accelerated upflow 
l)Caligari et al .1 fl995h . After emergence, the near-surface 
parts of the sunspots, pores, and other flux concentra- 
tions are rapidly cooled by radiative losses, which gives 
rise to an inward propagating cooling front accompanied 
by a local downflow. This leads to a decrease of the gas 
pressure and a concomitant intensification of the magnetic 
field in the upper layers of the flux tube, while the upflow 
from below and the downflow from above increase the gas 
pressure below the first few Mm depth. After a few days, 
the gas pressure has increased sufficiently to approach the 
ambient pressure somewhere between 2 and 10 Mm depth. 
A thin flux tube (requiring lateral balance of total pres- 
sure) would then lose its equilibrium and explode. A more 
general magnetostatic equilibrium could possibly be main- 
tained for some time but, in any case, the magnetic field 
becomes strongly weakened. Once the field strength falls 
below the equipartition value with respect to the exter- 
nal convective motions, the field becomes largely passive 
in that height range and progressively fragments, so that 
the upper, magnetostatic part becomes dynamically dis- 
connected from its roots. 

3. A simple model for dynamical disconnection 

In order to evaluate the explanatory potential of the sce- 
nario sketched above and to provide an illustrative exam- 
ple, we study the quasi -static evolutio n of a sunspot model 
very similar to that of lDeinzerl 1 19651) . 

3.1. Model description 
Evolution 

To realistically follow the sequence of events from flux 
tube emergence to disconnection would require a fully dy- 
namic realistic numerical simulation, which is beyond the 
scope of this paper. We restrict ourselves to sequences of 
quasi-static models of a sunspot and its underlying flux 
tube, which evolve in reaction to surface cooling and to 
the increasing pressure resulting from an inflow through 
the lower boundary of the model. Such an approach is jus- 
tified as long as the upflow velocity is significantly smaller 
than the Alfven speed, so that the magnetic structure re- 
mains almost unaffected by the dynamical pressure of the 
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flow. Starting from an isentropic initial state, we proceed 
in constant intervals of time and, for each time step, de- 
termine the entropy change due to radiative and convec- 
tive energy transport. Next we integrate the gas pressure, 
taking into account the growing base pressure owing to 
the assumed upflow, apply the corresponding adiabatic 
correction to the temperature, and calculate the modi- 
fied magnetostatic equilibrium. Then we move to the next 
timcstep. In the following subsections we explain the var- 
ious ingredients of the model in more detail. 

Magnetostatic equilibrium 

We basically follow the approa ch of iDeinzerl ll 19651) 
and c onsider a similarity solution CSch luter fc Temesvarvl 
119581) for an axisymmetric magnetostatic configuration 
with hydrostatic equilibrium along the field lines and a 
mixin g-length model of the convection zone l|Kiefer et alJ 
2000) representing the (fixed) external stratification. 

The similarity model assumes an axisymmetric mag- 
netic structure with a self-similar profile of the magnetic 
field as a function of the radial coordinate, r, in cylin- 
dric coordinates. The two components of the (untwisted) 
magnetic field are written in the form 



B z {r,z) = f(()B (z) 



B r (r, 



2 J ^> dz 



(1) 
(2) 



where Bq(z) is the magnetic field along the symmetry axis 
and C = T-J Bq{z). The function /(C), which describes the 
radial profile of the vertical field component, can be freely 
chosen. We follow the usual practice and use a Gaussian: 
/(C) = ex P( — C 2 )- Inserting the ansatz given by Eqs. (JJJ 
and PJ into the equations for magnetostatic equilibrium, 



= dp | Bz ( dBr dBz 

dr 47r \ dz dr 

Q = dp B r /6B r dB z 

dz 4n \ dz dr 



pa 



(3) 



(4) 



and integrating the first of these over r (from to oo) for 
constant z, yields the ordinary differential equation 



solving two ordinary differential equations, namely hy- 
drostatic equilibrium along the central field line of the 
axisymmetric configuration and an equation relating the 
field strength on the axis to the difference between the 
gas pressure on the axis and the (external) pressure far 
away from the axis. Note that the thin flux tube approx- 
imation with B 2 = y A = 8tt(p c — pi) is recovered in the 
limit $ — > 0. The left-hand side of Eq. J^J represents the 
contribution of magnetic stress integrated from the axis of 
the field configuration to infinity, which is not considered 
in the thin flux tu be limit. A detailed derivatio n of Eq. 
© can be found in lSchhiter fc Temesvarvl (|l958|) . 

We solve Eq. through an iterative relaxation pro- 
cedure, 



Vn+l 



e 

Vn + — 

Vn 



2-7T dz z 



vlVn+1 + 87r(p c -ft) 



,(7) 



where y n {z) is the nth iteration of the solution. This pro- 
cedure avoids numerical stiffness problems introduced by 
the term y 4 and allows us to specify boundary conditions 
for y on both sides of the computational domain. The 
prefactor y~ 3 has been introduced in order to increase the 
convergence and stability of the algorithm, which we op- 
timize by the choice of the free parameter e. 

Energy transport 

We consider vertical radiative transfer in the grey diffu- 
sion approximation. Since the energy flux of a sunspot 
umbra cannot be provided by radiation alone, a reduced 
level of convec tive energy transport is required in addition 
ilDeinzerlll965l) . The phenomenon of umbral dots is often 
taken as a manifestation of small-scale convection in the 
strong umbral magnetic field. In the absence of a better 
theoretical model, we treat the magneto-convective energy 
transport using a mixing-length approach, taking into ac- 
count the inhibiting effect of the magnetic field through a 
reduction of the mixing-length parameter. Horizontal ra- 
diative energy exchange can be neglected in the case of 
sunspots, which are much larger than the photon mean 
free path of the order of 100 km at optical depth unity. 
The time evolution of 7] then follows from 



2n y dz 2 



'II 



y -8tt(p c 



Pi 



where y = y/ Bq (z) , $ is the total magnetic flux, p c (z) = 
p(oo, z) is the (fixed) external pressure, and pi — p(0, z) 
is the pressure on the axis. The latter is simply given by 
hydrostatic equilibrium along the axis, viz. 



dpi 

dz 



Qi9, 



(6) 



where Qi(z) is the density profile along the axis and 
g(z) is the (depth-dependent) gravitational acceleration. 
Consequently, the self-similarity assumption reduces the 
determination of the magnetohydrostatic force balance to 



(5) . c ?Ii - - — (F ,+F ) 



dt dz 

where the radiative energy flux is given by 



ad 



16aT? dTi 

iKRQi dz 



(8) 



(9) 



The Rosseland mean opacity, Kn(gi,Ti), is interpolated 
from a table. Since the value of k,r significantly drops in 
the surface layers, the radiative heat conductivity reaches 
large values, so that we use a semi-implicit treatment of 
the radiative energy flux in order to avoid a too severe 
time step constraint. 
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The convective energy flux a long the flux tube is de- 
termined following ISoruitl Jl974h . 

i^conv = -bV^ (jf) V2 (J^ V2 6iC P T? /2 (V - V') 3/2 (10) 

where V is the logarithmic temperature gradient with re- 
spect to the logarithm of pressure, \x is the average molec- 
ular weight, TZ the gas constant, and l/H p the (constant) 
mixing length parameter as fraction of the local pressure 
scale height. The logarithmic temperature gradient gradi- 
ent V' reflects radiative energy exchange of the convective 
parcels and follows from 

V = V ad - 2u 2 + 2u(V- V ad + u 2 ) 1/2 , (11) 
with 

u =j-(±y 2 12 ^ 3 fv\ 1/a . (i 2 ) 

fy/a \H p J CpQiKftHp \ g J 

Here V a d denotes the adiabatic logarithmic temperature 
gradient and a is Stefan's radiation constant. The values 
of the geometric parameters used are a = 1/8, b = 1/2, 
and / = 3/2. 

In the case of a very strong systematic vertical flow 
with a time scale comparable to the turnover time of the 
convective motions, the mixing-length approach may be 
invalidated and, in addition, the advective entropy trans- 
port by the large-scale flow may become relevant. In our 
simulations, such a situation occurs only during the initial 
cooling phase when a rather strong transient downflow de- 
velops in the first 1-2 Mm depth (see the subsection on the 
implied vertical velocity below) . For the later evolution of 
the configuration, the mixing length model appears to be 
adequate for the purposes of an illustrative model. 

Equation of state 

We include the partial ionization of H, He, and Hc + . The 
values of density, heat capacity, and adiabatic temperature 
gradient as functions of pressure and temperature are de- 
termined in the course of the calculation by interpolation 
in a pre-compiled table. 

Boundary conditions 

For the temperature, we have to specify boundary condi- 
tions at both ends of the integration domain. At the upper 
boundary (2 = 0, taken to be the level of Rosseland op- 
tical depth unity of the external stratification) we fix the 
temperature to a value (typically 3500 K) that is lower 
than the temperature expected at the level r = 2/3 in the 
sunspot. At the bottom, we keep the temperature gradient 
adiabatic. 

For solving Eqs. JSJ and 10 we specify a boundary 
condition for p\ at the bottom of the integration domain 
and values for y at the bottom and at the top, respectively. 
The magnetic field strength is fixed at the upper boundary 



(typically at a value of 2000 G), while the thin flux tube 
relation, B = ^/8ir(p c — p{), is assumed to hold at the 
bottom. 

We integrate the hydrostatic balance upward from the 
lower boundary, specifying the value of the base pressure 
for each time step according to the influx of mass due 
to the assumed upflow. The effect of lateral expansion or 
contraction of the configuration is taken into account in 
the calculation of the base pressure. We can alternatively 
specify a constant upflow velocity, vq, or a constant total 
mass flux through the lower boundary of our integration 
domain. 

In the thin flux tube approximation, the total mass 
within a flux tube containing the (constant) magnetic flux 
$ is given by 

m = $ / — dz, (13) 
Jo b 

where zq represents the lower boundary of the integration 
domain. The base pressure, pi(zo), is adjusted each time 
step such that the change in total mass, Am, reflects the 
inflow of material across the boundary, 

Am = A{z )g i {z )v At^^^^-v At , (14) 

where A(z) = $/B(z) is the cross-sectional area of the 
flux tube. These relations also hold more generally for the 
self-similar solutions considered here, since the profile of 
the magnetic field in radius is the same at each depth, z, 
except for a scaling factor. Therefore, we can define an 
'effective' cross section of the configuration at each height 
on the basis of the axial field as A e g (z) = $>/B(z), where 
A e g depends on the radial field profile considered. A e s is 
in general different from the geometric diameter of the flux 
configuration (for instance, in the case of a Gaussian pro- 
file it formally has infinite extent). However, since the total 
mass of the flux tube and the mass flux across the lower 
boundary scale in the same way with the cross-section, 
the difference between A e s and the geometric cross-section 
has no influence on the solution. 

A change of the pressure at the base, Pi(zo), alters 
the mass within the flux tube through a change of q\ 
and B (the latter is the dominant contribution since 
(3 = 8npi/B 2 3> 1 except for the uppermost layers), which 
can be written formally as m = m[pi(zo)]. The adjustment 
of the base pressure Api(zo) has to be determined such 
that 

m[ Pi (z ) + A Pi (z )} = m[ Pi (z )} + $%^v At (15) 

B{zo) 

holds. Since the relation between Api(zo) and Am is non- 
linear [mainly through Eq. J^J] , we use a Newton iteration 
to determine the adjustment of the base pressure such that 
Eq. I|15(l is satisfied. Note that each iteration step requires 
a full integration of the Eqs. © and 10. 

For reasons of consistency with the self-similar solution 
it is required that the flow velocity has no variation over 
the cross section of the flux tube. This ensures that the 
magnetic field profiles remain self similar at all times. 
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Implied vertical velocity 

The inflow at the bottom, the radiative and convective 
cooling, and the resulting changes of the magnetostatic 
balance lead to a net change as well as to a redistribution 
of mass within the magnetic structure. In our quasi-static 
model, we can determine the implied vertical velocity on 
the axis by considering the equation of continuity: 

9 Qi d ( Qi\ 

d-tB + d- Z { V B) = °- (16) 
Starting with v — at the upper boundary, integration 
along the tube axis leads to: 

v{z) = ^777 / ^dz . 17 

Q\{z)otJ Q B 

We determine the time derivative as a finite difference 
from two consecutive time steps of the hydrostatic solu- 
tions. The feedback of the dynamic pressure on the solu- 
tion can be neglected as long as QiV 2 /2 <C p c — Pi- Even 
for the maximum inflow velocities that we consider (of the 
order of 1000 m/s) , the dynamic pressure is more than two 
orders of magnitude smaller than the difference between 
the external and internal gas pressures. In other words, the 
inflow velocity is at least one order of magnitude smaller 
than the initial Alfven speed. In all cases that we have 
considered, the hydrostatic approximation continues to be 
valid throughout the whole simulation. 

We have not considered the advection effects from the 
implied vertical velocity in the energy equation. In the 
deeper, homentropic part of the flux tube, such effects are 
irrelevant for the (adiabatic) temperature structure, while 
the downflow in the upper part could be temporarily im- 
portant, but rapidly diminishes to less than 20m-s _1 after 
the first phase of radiative cooling. The advection of low 
entropy by the downflow probably enhances the intensifi- 
cation of the magnetic field in the upper part (in effect, 
this is the convective collapse mechanism) . A quantitative 
evaluation of these effects requires a fully dynamic simu- 
lation, which is beyond the scope of this paper. 

Initial condition 

We initialize the problem with an isentropic stratification, 
using an entropy value equal to the entropy of external 
stratification at the lower boundary. Since the external 
stratification is strongly superadiabatic in the uppermost 
layers of the convection zone, the magnetic field strength 
at the lower boundary has to be chosen fairly large (around 
100 kG for a bottom depth of 12.5 Mm) in order to to 
avoid the loss of lateral pressure balance (i.e., an explo- 
sion) already in the initial state. This is a consequence of 
the magnetostatic approach chosen here. Therefore, the 
time evolution in the deeper layers as calculated here does 
not fully represent the real solar situation. On the other 
hand, the final state of an (nearly) exploded configuration 
is a fully consistent and relevant magneto-hydrostatic so- 
lution, which shows that a static sunspot is possible after 
dynamical disconnection. 




I , , , , I , , , , I , , , , I , , , ,_ 

0.0 0.5 1.0 1.5 2.0 

Depth (Mm) 

Fig. 1. Sequence of temperature profiles as a function of 
depth in the upper part of the model, showing the rapid 
inward propagation of a cooling front and the establish- 
ment of a largely stationary profile after a few hours. 
Zero depth corresponds to optical depth unity in the 
external stratification. The uppermost (thick full) curve 
shows the (time-independent) external temperature. The 
other curves show the temperature on the sunspot axis 
after 0.5h (dotted curve), lh (dashed curve), 4h (dash- 
dotted curve), 15 h (dash-triple-dotted curve), and 30 h 
(full curve), respectively, from the onset of radiative cool- 
ing (corresponding to the emergence of the flux tube at 
the surface). 



3.2. Numerical results 
Example run 

As an illustration of a typical evolution in the framework 
of our quasi-static model, Figs. ^CD show snapshots of 
temperature, magnetic field strength, and convective en- 
ergy flux, respectively, along the axis of the flux tube. 
In this case we have assumed a total magnetic flux of 
10 21 Mx and a constant mass flux at the bottom (located 
at 12.5 Mm depth) corresponding to an initial inflow ve- 
locity of 700 m-s -1 . Simulations of thin flux tubes give ris- 
ing sp eeds of that order at the same depth l|Caligari et alJ 
1995). The spacing of the numerical grid is 25 km. The 
temperature at the upper boundary (at optical depth 
unity of the external stratification) has been fixed at a 
value of 3500 K and the magnetic field strength at 2000 G. 
The ratio a = l/H p of the mixing length, I, to the lo- 
cal pressure scale height, H p , determines the total energy 
flux in the flux tube, which is rather well constrained by 
sunspot observations at 20-25% of the undisturbed so- 
lar flux. A value of a — 0.3 leads to an energy flux of 
about 22% of the undisturbed solar value at the top of the 
model. About 40 hours after the start of radiative cooling 
and bottom inflow, the field strength has fallen below the 
equipartition value (with respect to the convective veloc- 
ity in the exterior) at at depth of 4.7 Mm, implying the 
dynamical disconnection of the upper part from its roots. 
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Fig. 2. Snapshots from the time evolution of the magnetic 
field strength along the tube axis as a function of depth 
for the same case as shown in Fig. ^ (note the different 
depth range). The curves correspond to about lh (upper 
full curve), 10 h (dotted curve), 20 h (short-dashed curve), 
30 h (dash-dotted curve), 35 h (dash-triple-dotted curve), 
and 40 h (long-dashed curve) after the start of radiative 
cooling (emergence). The lower, thick full curve gives the 
equipartition field strength with respect to the convective 
velocities fro m the mixing lengt h model of the external 
stratification l|Kiefer et alJl2000|) . After an initial general 
drop, the field strength increases again down to a depth of 
about 4 Mm while a local minimum of the field strength 
develops and moves downward until, after about 40 hours, 
the field strength falls below the local equipartition value 
at at depth of 4.7 Mm. 
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Fig. 3. Profiles of the convective energy flux density (nor- 
malized by the undisturbed solar value) along the sunspot 
axis for the same instants of time as shown in Fig. |21 
(note the different depth range). Above the Wilson de- 
pression, the energy transport is taken over by the ra- 
diative flux (not shown here). The convective transport 
affects progressively deeper layers in order to supply the 
(almost constant) surface energy flux density of about 22% 
of the undisturbed solar value. The corresponding down- 
ward extension of the superadiabatic stratification leads 
to a growth of the region that is stabilized against dis- 
connection through the resulting increase of the magnetic 
field (cf. Fig. H 



At the same time, the inflow velocity has decreased to less 
than 50m-s _1 owing to the strong weakening of the field 
and the assumed constant mass flux. 

Fig. ^ shows the fixed external temperature stratifica- 
tion (thick full line) and five profiles of the temperature 
in the upper part of the model at various instants of time. 
After the emergence of the flux tube at the photosphere 
(start of the calculation), radiative cooling leads to the 
development of a cooling front, which rapidly progresses 
inward. Already after about 30 minutes, the level of opti- 
cal depth unity in the flux tube (the Wilson depression) 
has reached a depth of about 375 km and after a few hours 
the inward propagation of the profile has become very slow 
(dash-dotted curve, after 4 hours). The last curve (thin 
full line) shows the largely stationary temperature profile 
about 30 hours after the start of radiative cooling. The 
Wilson depression has reached a value of about 550 km 
with a temperature of about 4100 K at (Rosseland) opti- 
cal depth unity. 

For the same case, Fig. [5] shows snapshots of the mag- 
netic field profiles along the tube axis. The full curve shows 
the situation after about 1 hour with a steep inward rise 
of the field strength that reflects the not quite realistic 
initial condition. As time progresses, the depth gradient 
of the magnetic field strength flattens and a local mini- 



mum develops, which becomes deeper and moves down- 
ward while, at the same time, there is a slow increase of 
the field strength in the upper layers. Both effects are con- 
nected to the development of a superadiabatic stratifica- 
tion in the sunspot model, which results from the growth 
of the region with convective energy transport. This is 
illustrated in Fig. |3 which shows depth profiles of the 
convective energy flux density (normalized by the undis- 
turbed solar value) for the same instants as the magnetic 
field profiles in Fig. (note the different depth scale). The 
superadiabatic stratification in the region with significant 
convective energy flux leads to a decrease of the pressure 
scale height and an associated increase of the field strength 
and thus prevents the field from becoming weak in the up- 
per layers. As the superadiabatic region grows downward, 
the position of the field strength minimum also is shifted 
downward until the field strength eventually falls below 
the local equipartition level at a depth of about 4.7 Mm. 
The field strength at the bottom of the model has then 
reached about 10 kG, which is consistent with the results 
of thin flux tube simulations. The observable layers above 
the Wilson depression (optical depth unity, where the field 
strength remains around 2500 G) are largely unaffected by 
this whole development. 
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Fig. 4. Disconnection time (upper panel) and disconnec- 
tion depth (lower panel) as functions of the initial upflow 
velocity at the lower boundary. The mass influx is kept 
constant in time in the individual simulation runs (indi- 
cated by asterisks). The connecting curves represent spline 
interpolations. 



Dependence on parameters 

The simulation run presented in the previous subsection 
yields values for the time scale and the depth of the dis- 
connection that appear to be reasonable in view of what is 
indicated by observations of solar magnetic regions. How 
do these results depend on the various parameters and 
assumptions entering the model? 

The most important parameter of the model is the in- 
flow velocity at the bottom, which determines the rate 
of pressure buildup at the base. For the case of constant 
mass flux through the bottom, Fig. Ogives the disconnec- 
tion time and depth (determined as time and depth for 
which the field strength first equals the external equipar- 
tition value) as functions of the initial inflow velocity, all 
other parameters being the same as those of the reference 
run discussed above. The plots show that that smaller val- 
ues of the mass flux (smaller initial inflow velocity) lead 
to later and deeper disconnection. In the alternative case 
with constant inflow velocity (as opposed to constant mass 
influx), even much smaller inflow speeds of the order of 
100 m-s^ 1 lead to disconnection within few days and in a 
few Mm depth. 
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Fig. 5. Disconnection depth vs. disconnection time for 
various model runs. Plus signs: runs with constant mass 
influx (the same runs as in Fig. 0J; asterisks: runs with 
constant inflow velocity in the range 75 - 400 m-s -1 ; di- 
amonds: runs with different values for the initial pressure 
at the top, changing the initial mass in the flux tube by 
about ±10%; triangles: runs with different values of the 
total magnetic flux (10 20 Mx and 3T0 21 Mx, respectively); 
crosses: runs for which the magnetic field was determined 
according to the pressure difference (approximation of thin 
flux tubes). The curve illustrates the downward progres- 
sion of the zone with convective energy transport (cf. 
Fig-IHJ); it shows, as a function of time, the depth at which 
the convective energy flux density equals 10% of the undis- 
turbed solar value, about half of the surface flux density 
of the model sunspot. Except for the largest disconnec- 
tion depths (which may be affected by the proximity of 
the lower boundary of the model at 12.5 Mm), the points 
closely follow the curve, indicating that the disconnection 
depth is mainly determined by the downward progression 
of the convective cooling region. 

For high initial inflow speeds, the disconnection depth 
and time become largely independent of speed and reach 
values around 9 hours and 2.3 Mm, respectively. This satu- 
ration results from the low entropy of the gas in the upper 
layers of the model due to the strong surface cooling, so 
that a strong field is always maintained in these layers. In 
fact, switching off the radiative losses invariably leads to a 
dramatic weakening of the magnetic field at the very top 
of the integration domain (depth zero). 

We have varied a number of other parameters as well 
as initial and boundary conditions in order to evaluate 
their effects on the disconnection depth (z<j) and discon- 
nection time (id)- We find that the results are rather in- 
sensitive to the choice of temperature and magnetic field 
at the top (z = 0). The same is true for the (constant) 
value of the entropy in the initial stratification. In most 
cases we have taken the initial entropy equal to that of the 
external medium at the bottom of the model (12.5 Mm 
depth). Since the plasma in the flux tube originates from 
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deeper layers zone, its entropy could be somewhat larger. 
However, since nearly all of the entropy drop through the 
convection zone takes place above 12.5 Mm, even taking 
the entropy value corresponding to the bottom of the con- 
vection zone has almost no effect on the results. 

However, other parameters and conditions have a sig- 
nificant effect on the disconnection. It turns out that Zd 
and td decrease 

- with decreasing depth of the lower boundary, 

- with decreasing total magnetic flux, 

- with increasing initial mass content, 

- and also by keeping the inflow velocity (in contrast to 
the mass influx) constant in time. 

Likewise, determining the magnetic field by replacing the 
similarity solution by the thin flux approximation [taking 
the limit $ — > in Eq. (J5J] leads to more shallow and 
earlier disconnection. 

The quantitative results are combined in Fig. which 
gives Zd as a function of td for a number of cases. Although 
the various simulation runs represent quite different con- 
ditions and assumptions, the results suggest a universal 
relationship between disconnection time and disconnec- 
tion depth. The curve drawn in the figure is not a fit but 
gives, as a function of time, the depth at which the con- 
vective energy flux density equals 10% of the undisturbed 
solar value, about half of the emergent surface flux density. 
This curve therefore reflects the downward progression of 
the region of convective cooling as a result of the radiative 
losses at the surface (see Fig. EJ. 

The relationship between the disconnection depth and 
the downward progressing convective 'cooling front' is in 
accordance with our interpretation of the disconnection 
depth in the case discussed in the previous subsection. 
Without surface cooling, the strongest weakening of the 
field would occur near to the top of the model. The cool- 
ing leads to the development of a superadiabatically strat- 
ified layer connected with a strengthening of the magnetic 
field in the upper part of the tube, which grows downward 
as the surface flux has to be supplied by the thermal en- 
ergy stored in the deeper layers. This prevents the field 
in the upper layers from being weakened by the pressure 
buildup due to the inflow from below. The disconnection 
then takes place somewhere in the flank of the downward 
progressing convective flux profiles shown in Fig. |21 lead- 
ing to the relation between Zd and td apparent from Fig. El 
In the absence of convective energy transport, the cool re- 
gion would only extend down to the Wilson depression 
and the sunspot would be disconnected and disrupted im- 
mediately below that level, making it difficult to imagine 
that it could survive as a coherent entity thereafter. 

This interpretation is further supported by runs with 
varied efficiency of convection, which affects the downward 
progression speed of the region of convective cooling. For 
the case shown in Figs. QHHl we have a = l/H p = 0.3, 
which gives a heat flux density ratio of 22%. Larger 
(smaller) values of a lead to a faster (slower) downward 
extension of the convectively cooled region and to deeper 
(shallower) explosion: for a = 0.4 we find an explosion 



depth of about 6 Mm and a flux of 30% of the undis- 
turbed solar flux density, while a tube model with a = 0.2 
explodes at a depth of 2.3 Mm and shows only 14% of the 
undisturbed solar energy flux density. In both cases, the 
corresponding explosion depths and times are consistent 
with the downward propagation speed of the cooling flank. 

4. Discussion 

The results from our illustrative model indicate that dy- 
namical disconnection of emerged magnetic flux is possi- 
ble at a depth of less than about 10 Mm and within a 
few days after emergence if upflow velocities as shown by 
simulations of thin flux tubes are considered (of the order 
of a km/s at 10 Mm depth). However, we have to keep 
in mind that the assumption of a quasi-static state in the 
model forces us to use a too strong initial magnetic field 
in the deeper layers, so that the time evolution as shown 
here may not correctly represent what is actually happen- 
ing below a sunspot. A more realistic time evolution of 
the disconnection process requires a full MHD simulation. 
On the other hand, the disconnected state represents a 
fully self-consistent solution with realistic surface proper- 
ties and a hydr ostatic stratifi cation as indicated by ob- 
servations (e.e.. iBeckerd fl977i) . Our results show that a 
disconnected magneto-static sunspot model can be con- 
sistently constructed. 

We have taken the equipartition field strength with 
respect to the convective velocities in the mixing-length 
model of the (non-magnetic) external medium as the limit 
for the onset of disconnection. This is only a rough indi- 
cator and the detailed disconnection process certainly is 
more complex. For instance, we have not considered the 
possible effect of fragmentation by the interchange (flut- 
ing) instability at the peripher y of a flux tube bou nded 
by a current sheet l)Parkerll975(l . lMever et alJ l)l977|) have 
shown that buoyancy may stabilize a sunspot only down to 
a few Mm depth below the solar surface. Deeper down the 
destabilizing effect of the curvature of the surface bound- 
ing the flux tube dominates. In our case, the weakening 
of the field by the upflow leads to a hourglass-like shape 
of the sunspot, so that buoyancy effects in the lower part 
actually promote the instability. This could lead to frag- 
mentation even before the equipartition limit has been 
reached. 

The requirement for disconnection of emerged flux 
does not only apply to sunspots but covers all magnetic 
flux in bipolar regions. We suggest that the smaller flux 
tubes which form the plage and network regions are dy- 
namically disconnected by the same process that affects 
the sunspots. In fact, using the thin flux tube approxi- 
mation we find a shallow disconnection depth (of about 2 
Mm) for such structures. This could be even smaller since 
a large part of the radiative energy loss of a small tube is 
provided by lateral radiative heating in the surface layers, 
so that the cooling of the subsurface layers proceeds more 
slowly than in sunspots, leading to a more shallow discon- 
nection. It is well conceivable that the weakening of the 
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field by the upflow in small flux fragments extends up to 
the surface, so that all initial connectivity is lost already 
during emergence or shortly thereafter. The formation of 
intense fields and the patterns of flux distribution in plages 
then results from the interacti on with the local velocity 
fields at the solar surface (e.g., IVogler et aLll2005h . 

In which way does the fragmented flux below the dis- 
connected surface flux evolve further? Since the flux frag- 
ments are passive with respect to convection, they will be 
deformed and stretched by the convective motions, prob- 
ably keeping the field strength at about the equipartion 
value. However, this must not lead to a dynamical 're- 
connection' of the emerged part with its magnetic roots. 
We conjecture that this could be prevented by progressing 
fragmentation into ever thinner filaments, owing to various 
instabilities an d dynamical processes fe.g.. lSchusslerll 984: 
The fragmentation proceeds until the re- 
sistive length scale is reached, for which flux freezing is no 
longer valid and effici ent exchange of mass and heat (by 
radiation) is provided l|Schusslerlll987j) . The magnetic flux 
then has become truly passive and is kinematically trans- 
ported by the motions of the plasma, providing no fur- 
ther possibility for restoring the dynamical link between 
the emerged flux and its the deep-lying parent flux tube. 
Such a transition to passive, weak field after the explosion 
of a flux t ube has also been found in the numerical sim- 
ulations o f lRerrmel & Scniisslerl fcOOlh . Sustained outflow 
of buoyant plasma from the 'stumps' of the disconnected 
parent flux tube and the horizontal motion driven by the 
curvature force until a narrow U-loop remains could also 
provide a mechanism for the removal of the parent flux 
tube from its storage region at the bottom of the con- 
vection zone. These processes can actually be observed in 
thin-tube simulations (Caligari, Rempel & Schiissler, in 
preparation). 

5. Conclusion 

We have shown that, in the framework of a simple mag- 
netostatic model, the dynamical disconnection of emerged 
magnetic flux from its magnetic roots due an upflow along 
the field is a rather robust feature, which occurs over a 
wide range of inflow velocities and within a reasonable in- 
terval of time after emergence. The disconnection depth 
lies between 2 and 6 Mm for disconnection times less 
than 3 days, which is consistent with the observation that, 
shortly after the flux emergence phase, active regions do 
no longer reflect the dynamics of their magnetic 'roots' 
deep in the convection zone. The disconnection depth is 
mainly determined by the extension of the radiatively 
and convectively cooled subsurface layer at the time of 
disconnection. Radiative cooling and downflow lead to a 
strengthening of the field in the upper layers and thus pre- 
vent a loss of equilibrium at the surface, which would lead 
to a complete destruction of the visible sunspot. Our re- 
sults are consistent with findings from local helioseismol- 
ogy with regard to the subsurface structure of sunspots 



and a ctive regions fe.g..lZhao et al.ll20mUCouvidat et~ai] 
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